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ABSTRACT 

Diskoseismology, the theoretical study of normal mode oscillations in geomet¬ 
rically thin, optically thick accretion disks, is a strong candidate to explain some 
QPOs in the power spectra of many black hole X-ray binary systems. The exis¬ 
tence of g-modes, presumably the most robust and visible of the modes, depends 
on general relativistic gravitational trapping in the hottest part of the disk. As 
the existence of the required cavity in the presence of magnetic helds has been 
put into doubt by theoretical calculations, we will explore in greater generality 
what the inclusion of magnetic helds has to say on the existence of g-modes. 

We use an analytical perturbative approach on the equations of MHD to assess 
the impact of such effects. Our main conclusion is that there appears to be no 
compelling reason to discard g-modes. In particular, the inclusion of a non-zero 
radial component of the magnetic held enables a broader scenario for cavity non¬ 
destruction, especially taking into account recent simulations’ saturation values 
for the magnetic held. 
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1. INTRODUCTION 

There currently exists a rich structure in the power spectra observations of black hole 
X-ray binary systems, which includes high frequency (40-450 Hz) quasi-periodic oscillations 
(QPO). Relativistic diskoseismology, the formalism of normal-mode oscillations of geomet¬ 
rically thin, optically thick accretion disks, is a strong candidate to explain at least some of 
these QPOs. (For a review, see Wagoner 2008.) 

Diskoseismology’s perturbative approach assumes that the effects of magnetic fields have 
been incorporated in the background equilibrium solution and works with fluid perturbations 
in which magnetic helds play no effective role. The objective of this paper is to study 
analytically the effects of including small magnetic helds on the oscillations described by 
relativistic diskoseismology. 

Building on previous work (Fu & Lai 2009, hereafter FL), we use a local (WKB) analysis 
of the full MHD equations to examine how the magnetic held ahects the physics of radial 
wave propagation. The main diherence with FL is that we include all three components 
of the magnetic held, not just the vertical and toroidal cases separately. We do assume, 
however, that the toroidal magnetic held component is larger than the other components 
(using cylindrical coordinates r, 0, z). This assumption is supported by simulations (see 
Table 1). 

We show that diskoseismic g-modes are more resistant to magnetic-held disruption than 
previously thought. 


2. RELATIVISTIC DISKOSEISMOLOGY AND GRAVITATIONAL 

TRAPPING 

Within diskoseismology, some observed high-frequency oscillations in the outgoing radi¬ 
ation of black hole X-ray binary systems such as GRO J1655-40 are due to normal modes of 
adiabatic hydrodynamic perturbations. These modes are the result of gravitational driving 
and pressure restoring forces in a geometrically thin, optically thick accretion disc in the 
weak thermal, steep power-law state (Remillard & McClintock 2006). 

The study of diskoseismology reveals the existence of diherent types of oscillation modes. 
Of these, the fundamental g-mode (an axisymmetric inertial-gravitational mode that oscil¬ 
lates mainly in the vertical plane) is the strongest candidate for explaining one of the QPOs, 
being the most robust and observable: it lies in the hottest part of the disk, has the largest 
photosphere, and is located away from the uncertain physics of the inner boundary (Perez 
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et al. 1997). The p-modes, on the other hand, are less observable. They are only weakly 
affected by the magnetic fields (FL). 

This interpretation is not only supported observationally by peaks in the power spectral 
density, but the g-mode has been observed in hydrodynamic simulations as well (Reynolds 
& Miller 2009; O’Neill, Reynolds, & Miller 2009). 

As one can see in Fig. 1, the fundamental g-mode is trapped just under the maximum 
value of the radial epicyclic frequency K{r) in the absence of magnetic fields. Therefore, the 
existence of g-modes would be compromised by physical conditions that have an effect on 
the corresponding trapping curve. (The explanation of the other curves in the figure can be 
found in Section 4.) 


3. THE EFFECTS OF INCLUDING MAGNETIC FIELDS 

The inclusion of magnetic fields modifies the shape of the trapping curve in Fig. 1 and 
therefore threatens the existence of g-modes, given that it is not certain whether the inner 
boundary (the innermost stable circular orbit) can effectively trap this type of mode. 


3.1. MHD Equations 


The Newtonian ideal MHD equations for non-self-gravitating accretion disks are: 


| + V.(pv) = 0. 


( 1 ) 


with 


(9v 1 1 

— + (v ■ V)v = — vn - + -T, 

at p p 


dB 


V X (v X B), 


( 2 ) 

( 3 ) 


T^^(B-V)B. (4) 

These are seven equations for seven unknowns: the density p, the velocity v, and the magnetic 
field B. A barotropic pressure P = P{p) is assumed, <F is the pseudo-Newtonian gravitational 
potential, and of course 


V-B = 0. 


( 5 ) 
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We restrict our analysis to standard thin disks. We include the most important effects 
of general relativity by using the exact expressions for the orbital angular velocity f2(r) and 
the radial epicyclic frequency K{r): 

f2(r) = + a)“^, (6) 

K,{r) = f2(r)(l — 6/r + 8a/r^'^^ — . (7) 

There is then no need to specify <h(r, z). 

The standard approach includes the assumptions that the unperturbed background flow 
be axisymmetric, with v = rfl{r)(j), and that B = + B^z, where B^ is a constant. 

These forms of v and B satisfy the stationary {d/dt = 0) equations ([I]) and ([3]). The radial 
force balance is dominated by the centrifugal and the gravitational terms, while the vertical 
balance is dominated by the vertical pressure gradient and the gravitational term. In other 
words, the magnetic fields are non-dominant. 


3.2. Inclusion of Radial Magnetic Fields 

In order to extend the usual approach, we include small (compared to B^) radial mag¬ 
netic fields. The traditional exclusion of radial magnetic fields in analytical treatments has 
no justification other than aesthetic prejudice (since they destroy the solutions’ stationarity). 
In fact, simulations consistently yield radial fields which are larger than vertical fields (see 
Table 1). 

We now discuss what happens in the formalism once one introduces a radial B^r term 
in the magnetic field of an otherwise stationary system, in particular the one described in 
the previous subsection. If this term has the form Br = C/r, where C is a constant, then 
([5]) is immediately satisfied, while (jS]) yields: 

B;^-= Bf^ + AB4r,t), (8) 

where 

AB^{r,t) = rBr—t. (9) 

We note immediately that AB (r, t) B^ B^ for timescales of a few fluid oscillations 

(such that ilt ~ 1). 

The new B^ field satisfies equation ([2D, and (being parallel to v) yields 0 on the right- 
hand side of equation ((3D. 
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In the context of the present paper on perturbations, we take Vr as negligible. However, 
it is important to verify that such an assumption does not lead to inconsistencies at the un¬ 
perturbed (equilibrium) level. The important point is that once one accepts a non-vanishing 
Br, then one must allow for a non-zero Vr in order to preserve the equality on equation (P|); 
in particular, this equation acquires a 0-component. 

Fortunately, the implied value for v^. ~ rQ is in fact no larger than the one expected 
from standard (Novikov &: Thorne 1973) viscosity considerations alone, which is Vr ~ 

Here, a* is the “viscosity parameter” (shear stress/pressure), 

h disk thickness 

typical value of r ’ ^ ^ 

while br is defined in equation (IT^ . 


3.3. Perturbations 

We consider perturbations of Sp, Sv, and 5B oc Recall that our unperturbed 

magnetic field has the form 


B = B^{r)(j) + {C/r)f + B^z . 

We work with the assumptions 

Bl,Bl Bl<^ AttP ~ 47rpc^ 

(where Cg refers to the speed of sound) and define the following small parameter: 

h = = 

y/All pci Cs 


( 11 ) 


( 12 ) 


(13) 


where is the Alfven velocity. The linearized equations for the perturbations then become: 

- iu5p + --^{prSvr) + + -^{p6v^) = 0 , 

r or r oz 


c- c- ^ 1 ^ c-TT 1 

— tOJOVr — 2\ldVft) = (jr -“I- 

p par Anp 


—+ + ^ 5B^- 

r oz Or Or 


e- ^ V cT-r 1 

-tudvs+^dvr = Gs -<^ 11+1 — 

p pr Airp 


im d f d 1 

- B(j) -|- Bz-^—\-Br ( ^—I— 

r Oz \ Or r 



(14) 


27ipr 



(15) 

( 

dB, 


Airp \ r dr 

(16) 
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^ 1 ^ fTT 1 

— ioj5vz = Gz -7;-<^n H- 

p p dz Airp 


I R ^ I R ^ 

--D</) + -D^T^- r^rTT 

r oz or 


5Bz, 


iu}6Br = 
d 


(imBj 


d 


1 “t“ Bz r-. 1 5 Vr 

\ r oz 


imBr 


Sva 


Br SVz , 


d . „ a 

V(f) — 

imB 


d 


dfl 


iujSB^ = -^{B^6vr) + Bz-^Sv^ - B^—6vz+-^{Br6v^) + r—SB^ , 


— iCj6Bz = - - — (rSvr) 

r or 7 

We have used the dehnitions uj = uj — mfl and 


+ 


(■^^bS)Sv.. 

\ r Or 


G = ivn- il. 

p p 


( 17 ) 

(18) 

(19) 

( 20 ) 

( 21 ) 


3.4. WKB Analysis of Axisymmetric Oscillations 

We now restrict ourselves to WKB conditions, in which by dehnition all perturbations 
are cx Furthermore, we study for simplicity axisymmetric (i.e., m = 0) oscillations, 

which are also more observationally relevant. We use 511 instead of 5p. Assuming kr,kz S> 
1/r and B^ cx and using the dehnition p = dlnQ/d\nr, the following equations obtain: 

— -^5n + ikrSvr + ikzdvz H- ^^^^^(BsSB^+BzSBz+BrSBr) = 0 , (22) 

pci Anpci 

tk 1 Ji 

-^5n + iudvr + 2 Vl5vs + — (ikzBz+ikrBr) 5Br - —5Bs = 0 , (23) 

p A-np znpr 

iu6v^ - ^Svr + (ikzBz+ikrBr) 6B^ + = 0 , (24) 

212 dvrp dvrpr 

ik 1 

-^5n + iuSvz H- (ikzBz+ikrBr) 6Bz = 0 , (25) 

p dvrp 

ikzBzSvr + iu5Br—ikzBr5vz = 0 , (26) 

ikrB^dVr — {ikzBz+ikrBr) 5v^ + ikzB^dVz — iooSB^ — pQ6Br = 0 , (27) 

ikrBz5vr—ikrBr5vz — ioo6Bz = 0 . (28) 

The quantities Gz, G^, and Gr have been neglected. Being an odd function of z, Gz is 
negligible near the midplane z = 0, and goes away when vertically averaging. The fact 
that we are close to a purely axisymmetric conhguration (as discussed above) implies that 
G(f) is also negligible for our purposes, whereas the term containing Gr in equation flTSl) is 
smaller than the other ones because radial-force balance is dominated by the centrifugal and 
gravitational terms in thin disks. 
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4. CAVITY BEHAVIOR 


In order to study the behavior of the g-mode trapping cavity under the inclusion of 
magnetic helds, one hrst obtains a dispersion relation from the equations for the perturba¬ 
tions derived in the previous subsection. Once the characteristic equation has been obtained, 
one needs to isolate the appropriate branch for ca, which is the one that has k as its leading 
term when the magnetic held goes to zero. An exploratory way of doing this is working 
to zeroth-order in k‘^ (i.e., setting = 0). The curves in Fig. 1 were obtained with this 
assumption, by means of a numerical approach (using the values for bi from Table 1). From 
now on we work with a non-zero 


Before doing that, though, a word on dispersion relation branches. The branches that 
describe Alfven waves and slow magnetosonic waves, which are the ones responsible for the 
magneto-rotational instability (MRI), are different ones from the one we study in this pa¬ 
per (Balbus & Hawley 1998). In particular, the MRI branches are characterized by low 
frequencies and growth rates of order bz H, which can be smaller than the g-mode frequen¬ 
cies. More explicitly, the typical timescale for MRI growth tmri is related to the g-mode 
oscillation period Tg by the following formula: 


^ _ K 1/2 

n ^ 


(29) 


which is ~ 3 for typical values of the involved quantities. Even though there are MRI effects 
of shorter timescales (~ 1/H), these only occur at very short length scales (<C disk thickness). 
Furthermore, the (alpha model) viscosity induced g-mode growth timescale Tvisc (Ortega & 
Wagoner 2000) is related to tmri by: 


t'mri 

^visc 


a* 


(30) 


which means that the MRI might grow no faster than the viscous g-mode growth. These 
results indicate that the g-modes may survive in the presence of MRI driven turbulent eddies. 


Another potential source of g-mode disruption is given by energy “pumping” from short 
to long length scales observed in freely decaying MHD turbulence, “inverse-cascade” simula¬ 
tions (Zrake 2014). On closer inspection, however, it is reassuring to see that the oscillation 
frequencies produced by this mechanism at the g-mode length scales are in reality much 
lower than the g-mode frequencies. 


We now employ a perturbative approach in order to solve the problem. Recall that we 
work with the small quantities £ 1 and bl,b‘^ b'^ ^ 1, and that we assume m = 0. 

The perturbed MHD equations fl2^ -fl28 |) , to order bibj, lead to the relevant dispersion 





relation: 


where 


oj — ojq + Ahfjjhz + Xb(j)bj^ + , 


[Re(cno)]^ 


Re(A) 

Re(A) 

Re(r) 


Re ( 7 ) 
Re(/?) 




K^klclil + bl) 


{kl + + bl) - 




fcfcf + 2(1 — p + q)Q?- 
2VtrK^(k‘lcl — K^) 

{l-p + q)clnkr 3/2 ( ^ 


~ E 


n 


K 


TK^ 


pkl^klcl 




pVt^k, 


r 


2Q'^krkzCl 


~ e 

~' r) 


K 


Q 


K 


Q 


n. 

K 


Q. 


The leading imaginary contribution comes from the term 


Im(A) 


hKcll2(2 + p)Sl^ - K^] 
ia{ky, - k2) 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


the effects of which are small compared to the real-part terms. (We note that S> k‘^ by 
one order of magnitude.) 


We note that the implied inverse timescale for possible mode growth due to equation 
fl38|) is 1/r ~ b(j,bzE^^^il, which is much smaller than the one corresponding to purely viscous 
effects (no magnetic fields) on an fundamental g-mode, I/tvIsc ~ except for very small 
values of a*. 


We also note that the sign of Im(A) is not determined by our formalism as kr and kz 
could have either values of the sign. 


5. DISCUSSION 

We are now in a position to offer an improved assessment of the effects on diskoseismol- 
ogy of hnite magnetic fields, including the important radial magnetic helds. 

Our results can be best appreciated in a plot of the form shown in Fig. 2, which describes 
the behavior of the trapping cavity in terms of the vertical and radial magnetic helds. The 
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cavity is only destroyed (i.e., there is no value of the radius at which du{r)/dr vanishes) 
outside the corresponding ellipse, for sufficiently large and helds. This hgure was 
obtained by scanning the behavior of u{r) for different values of the Bi and determining 
where the cavity disappears. (The mode lives in the range of radii where > 0 for a given 
eigenfrequency.) 

In order to generate these results, the following ansatz was used: = ek"^, which is 

consistent with a radial mode size VhR [cf. eq. (5.1) in Perez et ah 1997], with e = 0.1, 
and with r] = 1 (as in FL), where 

f^±(^") = fl(r)(l — (39) 

is the vertical epicyclic frequency. In addition, we also used an alternative ansatz given by 
= e'^k‘^, corresponding to a larger radial mode size ~ R, motivated by the fact that the 
g-mode radial extension might increase as the concavity of u decreases. (This second ansatz 
gives the maximum radial g-mode extension that does not contradict the WKB assumption.) 

Note that the dependence on B^ is rather weak, as long as B^ 3> B^^Br- Importantly, 
within each ellipse, the maximum value of a;(r) does not typically change by more than 
about 15% (see Fig. 3 for a typical case), which means that the results are consistent with 
a constant QPO frequency within the present limits of observation. We should point out, 
however, that the numerical results of Fig. 1 imply a somewhat greater range of variation 
for the maximum value of ca, in potential disagreement with observations. (Recall, though, 
that these results assume k^ = 0.) 

Even though the perturbative results cannot be directly compared to FL (who study 
only the B = B^{r)(l) and B = BzZ special cases, separately), their results are consistent 
with ours in general terms. 

Our main conclusions in the present exploratory approach are the following. First, 
from the above discussion there seems to be no compelling reason to discard axisymmetric 
g-mode K trapping. While it is still true that the inclusion of magnetic helds modihes the 
cavity, the situation is not as devastating as implied by FL. Note in particular that the 
inclusion of a non-zero B^ potentially allows for slightly larger values of cavity-preserving 
\Bz\. Furthermore, the numerical results of Fig. 1 hint that the perturbative results may be 
underestimates of these \Bz\ values. 

More importantly, most simulations appear to produce values of B^ and Bz which he 
within or near each ellipse of Fig. 2. See Table 1 and corresponding bullets in Fig. 2. (Note, 
however, that there is an outlier point, not plotted.) 

In the second place, it must be kept in mind that possible diskoseismic explanations 
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of QPOs require only that the magnetic held be inside the ellipses in Fig. 2 during some, 
possibly small, fraction of the time, as the corresponding QPO duty cycles are observed to be 
much smaller than 100% (Remillard & McClintock 2006; Belloni, Sanna, & Mendez 2012). 


This work was supported by grant 829-A3-078 of Universidad de Costa Rica’s Vicerrec- 
toria de Investigacion and by grant FI-0204-2012 of MICITT and CONICIT. Travel funds 
provided by Stanford and Universidad de Costa Rica. 
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Table 1. 


Saturation Values of the Magnetic Field. 


Type of 3D Simulation 

Initial B Field 


\br\ 

h\ 

Reference 

shearing box, purely radial gravity 

vertical 

0.09 

0.15 

0.26 

1 

shearing box, purely radial gravity 

toroidal 

0.04 

0.07 

0.21 

1 

shearing box, purely radial gravity 

{B>=0 

<0.03 

<0.05 

<0.16 

2 

shearing box, including vertical gravity 

twisted toroidal 

0.09 

0.11 

0.28 

3 

shearing box, including vertical gravity 

vertical 

0.06 

0.08 

0.23 

4 

global, magnetically choked, H/R k. 0.5 

vertical 

0.08 

0.56 

0.56 

5 

global, H/R = 0.1 

toroidal 

<0.05 

<0.07 

<0.17 

6 

global, different temperature prohles 

vertical, weak 

<0.04 

<0.06 

<0.17 

7 


Note. — Saturation values of the magnitude of the magnetic held components are quoted, in 
the form of hi = VAi/cgi according to various 3D simulations. HjR refers to the the disk’s typical 
thickness to radius ratio. 


References. — (1) Hawley et al. 1995; (2) Hawley et al. 1996; (3) Shi et al. 2010; (4) Simon et 
al. 2011; (5) McKinney et al. 2012; (6) Parkin 2013; (7) Suzuki & Inutsuka 2014. 
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a=0 



r [GM/c^] 


a=0.95 



Fig. 1.— g-mode trapping cavity numerical estimates for different values of the magnetic 
field, assuming m = 0 and setting /c^ = 0 as a first approximation (finite kr results within a 
perturbative approach are presented in Fig. 2 and discussed in Section 5). The normalized 
magnetic helds bi are dehned in equation ([13]), and their values are taken from Table 1. (The 
listed bi values at the upper right corners correspond to the curves, in descending order.) 
For a given frequency, oscillating modes can exist below the corresponding curve. For high 
enough values of the magnetic held, the cavity is destroyed (not shown), i.e., the curve fails 
to have a maximum. Also shown is K{r), the leading term of the cavity in the absence of 
magnetic helds. Jagged curves are g-modes (shown here for the case of vanishing magnetic 
helds), dash-vertical lines mark the inner disk boundary at the ISCO. The upper and lower 
panels are for the respective cases of a = 0 and a = 0.95, where a = cJ/GM'^ is the black 
hole angular momentum parameter. (This hgure appears in color in the online version of 
the paper.) 
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Fig. 2.— Behavior of the g-mode trapping cavity as a function of the three (normalized) 
components of the magnetic held bi for different values of a. The cavity is preserved within 
each ellipse (dark and light for kj = k^e and k^ = respectively, corresponding to 

different radial mode sizes), and destroyed outside it. Extrapolations to perturbative analysis 
are indicated by dashes; they occur whenever the maximum of hi and hi is larger than 6^/4. 
Also shown are bullets corresponding to Table 1 simulation saturation values (or their upper 
bounds), but note that we plot ± the values, as they carry no sign. We do not plot the 
outlier point. (This hgure appears in color in the online version of the paper.) 
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Fig. 3.— Behavior of the variation of the maximum value of u within the (perturbative 
analysis) trapping-cavity ellipse for the case a = 0.8, 6^ = 0.1, kj = k^e in Fig. 2, as a 
function of bz and br- Different shades of gray represent percentual differences, from 0% to 
14%, with respect to the smallest maximum value. (See, however, the comment in the main 
text about the implications of the numerical results of Fig. 1.) 




